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ABSTRACT 

The first stars in the universe ionized the ambient primordial gas through various feedback processes. 
"Second-generation" primordial stars potentially form from this disturbed gas after its recombination. 
In this Letter, we study the late formation stage of such second-generation stars, where a large amount 
of gas accretes onto the protostar and the final stellar mass is determined when the accretion termi- 
nates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet 
(UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution cal- 
culations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates 
onto the star are ten times lower than in the case of the formation of the first stars. The lower 
accretion rates and envelope density result in the occurrence of an expanding bipolar Hii region at a 
lower protostellar mass M* ~ 10M Q , which blows out the circumstellar material, thereby quenching 
the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due 
to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates 
mass accretion onto the star at M* ~ 17 M . Although the derived masses of the second-generation 
primordial stars are systematically lower than those of the first generation, the difference is within a 
factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars 
are born as massive stars with tens of solar masses, regardless of their generations. 
Subject headings: cosmology: theory - early universe - stars: formation - stars: evolution - accretion 
- Hii regions 



1. INTRODUCTION 

Standard cosmology predicts that the cradles of the 
very first stars in the universe are massive (100 — 
1000 Mq) gas clumps formed in minih alos hundreds mil- 
lion years after th e Big Bang (e.g.. | Abel et al.l l2002t 
IBromm et all l2002t lYoshida et all 12003ft . Themassive 
primordial stars ionize their natal clouds by ultraviolet 
(UV) radiation. If these stars die without triggering su- 
pernova explosions, their relic Hn regions, where the ion- 
ized gas recombined subsequently, become a promising 
site for the formati on of second-generation stars (e.g., 
lOh fc Haimar][200l . The pre-stellar thermal evolution 
of these second-generation stars differs from that of the 
very first stars. The enhanced free electrons, catalysts 
for forming H2 molecules, in the relic Hii regions results 
in the higher H2 abundance than in the case of the first 
stars. The higher H2 abundance enables the production 
of cold (< 150 K) dense gas, where m ost of the deuterium 
is converted to molecular form (e.g.. lUehara fc Inutsukal 
20001 : [Nakamu ra fc Umemur a 2002; Naga kura fc Omukail 
20051 : Johnson fc Brommll2006f) by the exothermic reac- 
tion 



D+ + H 2 -> HD + H H 



(1) 
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With the additional radiative cooling via HD molecules, 
the temperature further falls to a few xlO K, which is 
close to the CMB temperature at high redshifts. By con- 
trast, in the case of first-star formation the temperature 
remains > 200K. Strong shocks or cosmic-rays can raise 
the ionization degree of a primor dial gas, thereby setting 



similar physical conditions (e . g.. iStacv fc Bromml 2007 
Safranek-S hrader e t al.l 120101 : llnayoshi fc Omukail 12011 : 
Omukail l2012f ) . The typical mass of the star-forming 
clouds is given by the Jeans mass when the contracting 
gas attains its minimum temperature 
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where the characteristic density and temperature are 
set by the HP cooling (e.g., lO'Shea et all 120051: 
lYoshida et all l2007al IMcGreer fc Brvanl 120081) . This 
cloud mass-scale is about an order of magnitude less mas- 
sive than that for the formation of the very first stars: 
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where the characte ristic density and te mperature are set 
by the H 2 cooling (iBromm et al.lll999() . 

Because of the difference in the mass-scale of their na- 
tal clouds, the first- and second-generation primordial 
stars have been classified into different sub-populations: 
Popu lation III. 1 (Pop II I.l) and III.2 (Pop III.2) stars 
fe.g.. lO'Shea et al.ll2008h . Until recently, the standard 
picture of primordial star formation postulated very mas- 
sive (a few xlOO M Q ) Pop III.l stars and less massive 
Pop III. 2 stars (a few x 10 M Q ). However, their final stel- 
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TABLE 1 
Included Deuterium chemistry 



No. 


Reactions 


References 


K 1 


F)+ 4- p \ F) -t- 'V 




K2 


D + H+ -> D+ + H 


1 


R3 


D+ + H -> D + H+ 


1 


R4 


D + H -> HD + 7 


2 


R5 


D + H 2 -)• H + HD 


3 


K6 


D+ + H 2 -> H+ + HD 


2 


K7 


HD + H -> H 2 + D 


2 


R8 


HD + H+ -> H 2 + D+ 


2 


R9 


HD + 7 -> H + D 


1 


Note. — 


1: Galli & Palla (1998), 2: 


Standi ct al. (1991), 3: 



lar masses are still not clear. Both in the Pop III. 1 and 
III. 2 cases, very small (~ 0.01 M Q ) embryo protostars 
form after the gravitational collapse of the nata l clouds 
(e.g.. IQmukai fc Nishi|[l99l lYoshida et al.ll2008D . These 
protostars then rapidly grow in mass by accretion of gas 
from their surrounding envelopes. The final stellar mass 
is thought to be set when the mass accretion ceases. 

The key physical mechanism that controls the final 
stellar mass is stel lar radiation feedback against the ac- 
cretio n flow (e.g.. jMcKee fc Tanl 120081 : IHosokawa et al.l 
1201 ll: IStacv et all 120121) ! In our previous work 
(|Hosokawa et al.ll2011l hereafter Hll), we have reported 
the first radiation-hydrodynamic simulations for the Pop 
III. 1 case that cover the long-term (~ 10 5 years) evo- 
lution after the birth of the embryo protostar. We have 
demonstrated that the strong stellar UV feedback indeed 
terminates protostellar accretion, and shown that the re- 
sulting final stellar mass is ~ 40 M®. This mass range of 
Pop III. 1 stars is much lower than previously postulated 
and is rather close to that expected for Pop III. 2 stars. 

In this Letter, we apply the same method to a typical 
Pop III. 2 star-forming cloud f ound in a cosmological sim- 
ulation (|Yoshida et al.ll2007bD . We find that although the 
final masses of Pop III. 2 stars are systematically lower 
than those of Pop III. 1 stars, the difference is by a fac- 
tor of only a few. This result suggests that, whereas 
lower-mass stars potentially form as satellite-like objects 
via g ravitational fragmentation of the star-forming cloud 
(e.g^ [Machida et al.N2008t iTurk et "aTH2009t IStacv et al.l 
120101: iClark et al.ll201l[) . the majority of primordial stars 
are tens of solar masses regardless of their generations. 
The typical stellar systems in the early universe are likely 
to be such massive stars accompanied by a small number 
lower-mass satellite stars. 

2. NUMERICAL METHOD 

We employ the numerical method developed in Hll 
with necessary extensions. The hydrostatic evolution 
of the central protostar and the evolution of the ac- 
creting primordial gas are consistently solved (see Hll 
for details). The dynamics of the accretion flow is 
calculated with a two-dimensional (2D) axisymmetric 
radiation-hydrodynamics code (e.g.. lYorke fc Sonnhalterl 
120021 lYorke fc Wel9 119961 ). In addition to the chemical 
reactions implemented in Hll, we also solve the deu- 
terium chemistry with three species of D, HD, and D + 
(Table 1). The HD-line cooling rate is calc ulated consid- 
ering the first four rotational levels as in iGalli fc Pallal 



(1998). The central protostar is replaced with a sink 
cell, whose size is ~ 10 AU. This masking of the 
very vicinity of the star, where the gas density exceeds 
10 11 cm~ 3 , enables us to follow the long-term evolution. 
The mass accretion rate onto the protostar is evaluated 
from the mass influx to the sink cell. With the accre- 
tion rate provided, the protostellar evolution is calcu- 
lated by solving the stellar interior structure numericall y 
(e.g.. IQmukai fc Palldl200H: IHosokawa fc Omukaill2009ft . 
The stellar evolution then gives the stellar luminosity 
and effective temperature to the radiation hydrodynam- 
ics, which follows the interplay between the accretion 
flow and stellar radiation. 

As in Hll we calculate the dynamics of the accre- 
tion flow within 0.3 pc of the protostar. The initial 
ga s distribution i s config ured with the simulation output 
of Yoshid a et al.l (|2007bl ). who study the formation of a 
protostellar core in a relic Hn region, performing three- 
dimensional (3D) cosmological simulations. We reduce 
their 3D data at the moment when the central density 
reaches ~ 10 14 cm" 3 to 2D axisymmetric data by averag- 
ing over azimuthal angles. The original 3D data has some 
non-axisymmetric structure, but the deviation from the 
axial average is not large within the computational do- 
main. Moreover, deviation from the axi-symmetric struc- 
ture is smaller in the inner part of the cl oud. Our com- 
putatio nal domain contains ~ 150 M©. lYoshida et al.l 
(|2007allbh show that in the Pop III. 2 case the mass of the 
cloud that experiences the dynamical run-away collapse 
is about 40 M Q . In our initial condition, this part of the 
cloud is contained within 0.05 pc of the protostar. In 
comparison with the Pop III. 1 case (Hll), the density 
in the accreting envelope is lower, reflecting the lower 
temperature in the pre-stellar collapse stage. The grid 
resolution and boundary conditions are the same as in 
the fiducial case studied in Hll. 



3. RESULTS 

Immediately after the onset of the calculation, the stel- 
lar mass increases rapidly by accretion via a circumstcl- 
lar disk. The density within the disk is high enough 
(> 10 8 cm -3 ) for the three-body reaction to convert 
most of the H atoms to H2 molecules. This molec- 
ular disk extends to ~ 500 AU from the star when 
the stellar mass reaches 5 Mq. The mass accretion 
rate onto the protostar decreases with increasing stel- 
lar mass, e.g., 1.7 x 10~ 3 M yr" 1 at M* = 2 M Q , while 
6.1 x 10~ 4 Mq yr _1 at M* = 5 M Q . These rates are 
about an order of magnitude lower than those in the 
Pop III.l case, due to the lower density in the accret- 
ing envelope resulting from the lower gas temperature 
during the run-away collapse stage: for the collapse of 
a self-gravitating cloud, the accretion rate is given by 
M» ~ c 3 /G oc T 15 , where c s and T are the s ound 
speed and gas temperature in the accreting envelope (jShul 
[19771) . In the Pop III.2 described above, the gas 

temperature is lower than in the Pop III.l case because 
of the efficient radiative cooling via extra H2 and HD 
molecules. 

Figure QJ(a) shows that, early on, the protostellar ra- 
dius is 30 — 40 Rq and decreases slightly as the accretion 
rate decreases. The protostellar evolution can be un- 
derstood by comparing the two timescales, namely the 
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Fig. 1. — Evolution of the protostellar radius (upper panel) and 
ionizing (hu > 13.6 cV) photon production rate (lower panel). The 
blue and red lines represent the fiducial Pop III.2 and III. 1 cases 
in each panel. The black dotted line in the upper panel shows 
the mass-radius relation of metalless zero-age main-sequence stars. 
The open and filled circles mark the characteristic epochs of the 
protostellar evolution, beginning of the KH contraction and the 
arrival of the protostar to the ZAMS. The filled diamonds in the 
lower panel denote the epochs when the Hu region begins to grow 
toward the polar directions of the disk. 



stellar Kelvin-Hclmholtz (KH) timescale 

GM*M* 



and the accretion timescale 



M ' 



(4) 
(5) 



The former is the times over which the star radiates away 
its internal heat content, whereas the latter is the time 
over which the stellar mass doubles by accretion. Ini- 
tially, the KH timescale is much longer than the accretion 
timescale, i.e., the star has no time to radiate its internal 
heat before its mass grows by accre tion. This is the so - 
called adiabatic accretion stage fe.g. JStahler et aH [l986). 
Since the stellar luminosity L* rises rapidly with increas- 
ing stellar mass, the KH timescale eventually falls below 
the accretion timescale at M* ~ 6 Mq, where the stellar 
radius marks its maximum. After this the protostar loses 
its internal heat efficiently by rad iative diffusion and con - 
tracts quasi-hydrostatically (e.g.. lOmukai fc Pa lla 2003). 
During this KH contraction stage the stellar effective 
temperature T e s varies as oc R~ 2 and H ionizing-photon 
(hv > 13.6 eV) luminosity therefore increases signifi- 
cantly (Fig. []}b). 

After the stellar mass exceeds 10 Mq, an Hu region 
begins to develop, growing in the polar directions of the 
disk. Figure [2} (a) shows the gas distribution in the 




Fig. 2. — The structure of the accreting envelope surrounding the 
protostar. The panels (a) and (b) show the snapshots when the 
stellar mass is 14.5 Mq and 16.5 Mq, respectively. The elapsed 
time since the birth of the embryo protostar is also indicated on 
the upper-left corner. Each panel shows the spatial distributions 
of temperature (upper left), density (upper right), and molecular 
fractions of hydrogen (lower left) and deuterium (lower right). The 
arrows represent the velocity field with different scalings between 
the upper and lower panels. The open triangles in the lower panel 
indicate directions only for velocities higher than 10.5 km/sec. The 
dashed lines in the lower half of each panel define the edge of the 
Hu region. 



vicinity of the protostar at M* = 14.5 Mq. We see 
that the bipolar Hil region is breaking out from the ac- 
creting envelope. H2 and HD molecules are completely 
photodissociated by stellar H2 dissociating (11.2 eV < 
hv < 13.6 eV) photons except in the shadow of the disk. 
Within the Hu region the ionized gas flows outward at a 
velocity of several 10 km/s. The high pressure in the Hu 
region drives a strong D-type shock, which compresses 
and heats up the surrounding materials to ~ 10 3 K. The 
shocked gas has an outward velocity of several km/s. The 
Hil region dynamically expands into the accretion en- 
velope and its opening angle increases with protostellar 
mass (Fig. U-b). 

The growth of stellar mass via accretion continues even 
after the formation and initial expansion of the Hu re- 
gion. However, the accretion rate is reduced by radiative 
feedback when compared to the no- feedback case (Fig.[3j . 
In contrast to the stellar UV radiation, the shock driven 
by the expanding Hu region propagates around the disk 
and reaches the accreting envelope in the disk's shadow. 
The shocked gas flowing outward hinders accretion from 
the envelope onto the disk. The isolated disk exposed 
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Fig. 3. — Evolution of the mass accretion rates onto the proto- 
star. The blue solid curve represents the fiducial Pop III. 2 case. 
The asterisks on this curve indicate the moments of the snapshots 
presented in Figure [2] The open and filled circles mark the char- 
acteristic epochs as in Figure [T] The blue dotted curve denotes 
the reference case, whereby the stellar UV feedback was switched 
off. The red solid and dotted curves represent the fiducial Pop 
III. 1 cases with and w ithout the stellar UV feedback taken from 
IHosokawa eTall poll ). 

to the stellar UV radiation thins down via photoevap- 
oration. Accretion from the disk onto the star finally 
terminates at M* ~ 16.6 M©. This coincides with the 
moment when the protostar reaches the zero-age main- 
sequence stage (Fig.QJb). 

Although the above overall evolution is qualitatively 
similar to that of Pop III. 1 stars (Hll), the resultant 
final stellar mass M* ~ 16.6 M© is lower than for the 
fiducial III. 1 case ~ 43 M© (Figure [3]). This is mostly 
due to the fact that the protostellar evolution differs 
with the different accretion histories between the Pop 
III. 1 and III. 2 cases (Fig. [T]). At the lower accretion 
rate, in general, the KH timescale becomes comparable 
to the accretion timescale at a lower luminosity. The 
protostar thus begins to lose its i nternal energy and con- 
tracts at lower stellar mass (e.g., IQmukai fc PaIial[2Q0l 
IHosokawa feO mukal 2009). The stellar ionizing lumi- 
nosity significantly increases during the KH contraction 
stage. Therefore, the UV feedback is effective at a lower 
stellar mass in the Pop III. 2 case. Although the mecha- 
nism of dynamical expansion of an Hll region and sub- 
sequent photoevaporation of the disk, which ultimately 
shuts off the mass accretion, is common between the Pop 
III. 1 and III. 2 cases, the required stellar ionizing lumi- 
nosity for this mechanism to operate is much lower in 
the Pop III. 2 case than for the Pop III. 1 case (Fig. [TJ- 
b). This is because the Hn region expands more easily 
through the lower-density Pop III. 2 accretion envelope. 
The accompanying protostellar outflow and photoevap- 
oration rates are also smaller than in the Pop III. 1 case. 
In the Pop III. 2 case, however, the weak UV stellar feed- 
back is enough to shut off the slow mass accretion toward 
the protostar. 

4. DISCUSSION 

In this Letter, we have shown that radiative feedback 
limits the final mass of a "typical" Pop III. 2 protostar 
to < 20 Mq. With all the possible variations of ini- 
tial conditions among different primordial gas clouds, we 
naturally expect that the resulting final stellar masses 
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Fig. 4. — Evolution of the mass accretion rates onto the protostar 
with different parameters. The green, blue, and black lines repre- 
sent the cases with different a- viscosity parameters of ao = 0.3, 0.6 
(fiducial case), and 1.0. The magenta lines denote the test cases 
where the rotation of the cloud is artificially reduced to 30% of the 
original values. The solid and dashed lines represent the evolution 
with and without the stellar UV feedback in each case. 



would have some distribution. Cosmological simulations 
predict that, in the Pop III. 1 case, the accretion rate onto 
a star could differ by more than two orders of magnitude 
among different star-forming clouds which form through 
vario us halo assembly histories (e.g., IQ'Shea fc Normanl 
120071 ). Although such systematic studies have not been 
conducted for the Pop III. 2 case, we expect a similar 
variance and hence some distribution of the final stellar 
masses. Here we examine potential variances of the final 
stellar masses in our simulations. 

For instance, we use the a-viscosity formalism 
(jShakura fc Sunyaevll 19731 ) for mimicking the torque ex- 
erted by non-axisymmetric sp iral arms formed in massive 
self-gravitating dis ks (e.g., lYorke fc Sonnhalterl 120021 : 
iKuiper et al.l l20T0). Recent 3D simulations show that 
the e ffective a-parameter is around 0.1 — 1 in this case 
(e.g. JClark et al.l[201 It IKuiper et al.ll2011[ ): we adopt an 
a-parameter in the mid-plane ao = 0.6 (see Hll) for 
the fiducial case discussed above. Figure [4] shows the 
accretion histories and resulting final stellar masses for 
different values of ao- Analogously to Hll, the final stel- 
lar mass increases with the value of ao: ~ 13.6 M© 
for a = 0.3, whereas ~ 20.0 M© for a = 1.0. Fig- 
ure |4] also shows the cases where the initial rotation of 
the natal cloud is reduced to 30% of the fiducial value. 
For this case the stellar UV feedback operates only af- 
ter 20 M Q of the gas has accreted onto the star and the 
resulting final stellar mass is M* ~ 28 M©. 

Although this mass range is systematically lower than 
that of Pop III. 1 stars derived in Hll, the difference 
is within a factor of a few. In reality, the existence 
of photodissociating background radiation could further 
erase the difference in the final masses of Pop III. 1 and 
III. 2 stars. Recall that the lower masses of Pop III. 2 
stars originates from the enhanced abundances of H2 
and HD, and the resultant more efficient cooling during 
the pre-stellar collapse of the natal cloud. However, the 
presence of even a weak photodissociating background 
(~ 10 _22 erg cm~ 2 sr _1 sec -1 Hz -1 ) inhibits the forma- 
tion of H2 and HD molecules, thereby making the ther- 
mal evolution almost the same in the Pop III. 1 case (e.g., 
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lYoshida et an2007bl : IWoIcott-Green fc Haimanll20TTI) . If 
this is the case, the mass scale of Pop III. 1 and III. 2 stars 
would be identical. 

Although not considered in this study, gravitational 
fragmentation of the disk would also contribute to shape 
the mass distribution of primordial stars. For the 
Pop III. 1 case, recent 3D numerical simulations show 
that the circumstellar disk easily becomes gravitation- 
ally unstabl e to fra gmentation ( e.g.JMachida et alj|2008t 
IStacv et all 120101: Iciark et all I201H IGreif et al.l 12012ft . 
However, most of these studies consider the early evolu- 
tion only, far before the mass accretion ceases and hnal 
stellar masses are determined. One can naively expect 
that disk fragmentation reduces the final stellar mass 
as the accreting material is shared by multiple proto- 
stars born from the fragm ents (so-called " fra gmentation- 
induced starvation", e.g. JPeters et al.ll2010h . 

The nature of mass accretion with disk fragmentation 
might be even more complex during later evolutionary 
phases. For instance, gravitational interactions between 
the multiple protostars could redistribute angular mo- 
mentum. If angular momentum is carried away by a 
small number of (proto-)stars thrown out of the system, 
accretion of the remaining gas onto the remaining star(s) 
could be enhanced. Moreover, merging of protostars in 
the disk is also possible. The net result is that mass ac- 
cretion onto the star would be strongly time-dependent 
or even have the consequence of stochastic burst-like ac- 
cretion events which potentia lly change the evo lution of 
the protostars (but see, e.g.. iSmith et al.ll2012h . These 
effects should be examined in future long-term 3D simu- 



lations. 

If some of the low-mass stars born via the disk frag- 
mentation survive without significant mass growth, they 
would be low-mass satellite-like stars orbiting around the 
massive stars whose final masses are limited by stellar 
radiative feedback. Future work would explain how this 
effect differs between the Pop III.l and III. 2 cas e s (e.g ., 
iMachida et ail I2009H . As shown in lClark et all (|201ll) . 
disk fragmentation is controlled by a balance between 
the mass transfer rates onto and through the disk. Sys- 
tematic difference in accretion rates onto the star-disk 
systems among the Pop III.l and III. 2 cases would result 
in basic differences in the effects of disk fragmentation. 

Our results provide a new picture that the typical mass 
of the primordial stars are tens of solar masses regardless 
of their generations. This explains why the signatures 
of pair-instability supernova, which is the final fate of 
140 M Q < M» < 260 Mq stars, have not been seen in the 
abundance patterns o f Galactic metal-po o r stars (e.g., 
Tumli nson etall [20p: iFrebel et al.l [20091: iCaffau et all 



201 11) . 
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